Journal  of  Power  Sources  184  (2008)  622-626 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


pri 

SllbiautS 


Short  communication 

Thermal  behavior  of  overcharged  nickel/metal  hydride 
batteries 

Dahe  Li,  Kai  Yang*,  Shi  Chen,  Feng  Wu 

School  of  Chemical  Engineering  &  Environment,  Beijing  Institute  of  Technology,  Beijing  100081,  China 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  16  January  2008 

Received  in  revised  form  15  March  2008 

Accepted  18  March  2008 

Available  online  29  March  2008 


Keywords: 

Ni/MH  battery 
Thermal  model 
Microcalorimeter 


This  work  provides  a  two-dimensional  thermal  model  for  cylinder  Ni/MH  battery.  Thermal  model  is  devel¬ 
oped  to  analyze  the  thermal  behavior  of  the  battery  when  charged  and  overcharged.  Quantity  of  heat  and 
heat  generation  rate  of  the  battery  during  charge  and  overcharge  period  are  studied  by  quartz  frequency 
microcalorimeter.  Heat  generation  curve  is  fitted  into  a  function,  and  heat  transport  equation  is  solved. 
Analysis  with  the  model  and  experiment  show  that  temperature  rise  is  about  3  °C  and  difference  between 
the  model  and  the  experiment  is  no  more  than  0.1  °C. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electric  vehicles  (EV)  and  hybrid  electric  vehicles  (HEV)  are 
becoming  more  attractive  because  of  growing  demand  for  power 
sources.  EV  and  HEV  require  the  battery  with  high  specific  energy, 
high-rate  capability,  long  cycle  life,  and  low  environmental  impact 
as  power  storage  device,  which  boosts  development  of  battery 
technology.  At  present,  most  research  concentrated  on  Ni/MH  and 
lithium/polymer  battery  because  the  two  kinds  of  batteries  are 
green  and  have  good  performance.  Their  performances  are  both 
sensitive  to  temperature  change,  temperature  rise  even  caused 
some  safety  problems  of  lithium/polymer.  Therefore,  Ni/MH  bat¬ 
tery  is  a  popular  choice  for  EV  and  HEV  power  storage  device. 
In  order  to  avoid  potential  problems  caused  by  temperature  rise, 
information  about  heat  generation  is  required. 

It  is  reasonable  to  develop  a  thermal  model  for  the  battery  to 
analyze  the  impact  on  performance  when  temperature  is  rising 
owing  to  the  necessity  of  a  thermal  management  for  the  battery 
[1-3].  Several  methods  are  used  to  simulate  the  battery  behavior 
like  computational  fluid  dynamics  (CFD)  and  finite  element  meth¬ 
ods  (FEM)  [4].  Bernadi  et  al.  [5]  have  set  up  a  general  energy  balance 
for  battery  systems  with  assumption  of  uniform  heat  generation. 
After  that  Chen  and  Evans  [6,7]  developed  several  two-dimensional 
and  three-dimensional  thermal  models,  and  presented  some  cal¬ 
culation  methods  for  model  parameters.  Shi  et  al.  [8]  studied  rapid 
charging  of  spirally  cylindrical  nickel/metal  hydride  battery  with 
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thermal  model.  Sato  et  al.  [9,10]  examined  many  heat  genera¬ 
tion  ways  for  electric  vehicle  battery  and  found  that  the  battery’s 
temperature  rise  was  usually  caused  by  heat  generation  due  to 
electrochemical  reactions  and  Joules  effect.  Wu  et  al.  [11]  vali¬ 
dated  the  cooling  effect  of  auxiliary  cooling  device  with  thermal 
model  and  showed  the  state  of  charge,  open-circuit  voltage,  inter¬ 
nal  resistance,  power  and  available  energy  changed  with  battery 
temperature  rise.  The  mathematical  simulation  of  heat  transport  in 
the  battery  has  high  efficiency  and  fairly  low  cost  compared  to  labo¬ 
ratory  experiments,  moreover,  the  agreement  between  theoretical 
and  experimental  values  was  good  in  many  earlier  works. 

Once  the  battery  is  overcharged,  oxygen  generated  at  positive 
electrode  reacts  with  hydrogen  at  negative  electrode,  and  consid¬ 
erable  amount  of  reaction  heat  accompanies  recombination.  The 
working  condition  for  the  battery  would  be  much  worse  if  the  bat¬ 
tery  were  overcharged.  As  volume  and  weight  in  a  vehicle  are  rather 
limited,  the  battery  system  has  to  be  smaller  in  order  to  take  up  less 
space.  Now  that  temperature  is  an  important  factor  of  the  battery 
performance,  the  specific  cooling  device  must  be  effective  to  avoid 
excessive  temperature  rise.  We  will  discuss  heat  generation  in  both 
normal  charge  and  extreme  charge  and  develop  a  model  for  the 
whole  process  and  discuss  the  model  and  the  experiment  results. 
We  hope  that  some  problems  can  be  found  like  abnormal  hot  spots 
or  areas  to  improve  thermal  management. 

2.  Experiments 

The  experiment  uses  cylindrical  Ni/MH  battery  with  spirally 
design,  rating  capacity  of  8  Ah  and  actual  capacity  of  7.5  Ah  after 
cycles.  The  battery  includes  electrodes,  electrolyte  and  separator 
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Table  1 

Some  parameters  of  tested  battery 
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Parameters 

Values 

R,  radius  (m) 

0.01609 

H,  height  (m) 

0.0605 

M,  mass  (kg) 

0.18909 

V,  volume  (m3) 

0.000049205 

and  is  assumed  to  be  axial  symmetric.  Some  details  are  in  Table  1. 

The  tested  battery  is  installed  in  a  special  device,  which  is 
designed  to  protect  the  battery  from  short  circuit.  During  the  exper¬ 
iment,  a  microcalorimeter  with  a  quartz  frequency  thermometer  is 
used  to  measure  the  surface  temperature  of  the  battery  and  the 
heat  generation  rate.  The  whole  set  is  in  the  microcalorimeter,  as 
shown  in  Fig.  1 .  All  data  are  recorded  by  a  PC. 

The  Ni/MH  battery  is  charged  to  fixed  SOC  (State  of  Charge)  in 
1 C  rate  first:  0, 30, 50, 70, 90,  and  100%.  Then  the  battery  is  charged 
in  1C,  3C,  and  5C  rate  to  150%  of  its  actual  capacity.  The  charging 
process  is  controlled  by  an  Arbin  instrument,  data  are  saved  by  the 
program  on  a  PC. 

3.  Results  and  discussions 

3.1.  Experiment  results 

The  quartz  frequency  thermometer  shows  the  surface  temper¬ 
ature  rises  after  1C,  3C,  and  5C  rate  charge  are  2.714,  2.883,  and 
2.826  °C.  Heat  generation  rates  are  shown  in  Figs.  2-4.  Curves  in 
the  three  figures  are  charged  in  1C,  3C  and  5C  rate  from  SOC  of 
0,  30,  50,  70,  90,  and  100%,  respectively.  It  seems  that  in  a  certain 
charge  rate,  a  curve  is  moving  to  the  left  as  SOC  increasing  from  0 
to  100%. 

During  normal  charge  process,  heat  generation  rate  is  slowly 
increasing.  Joule  heat  plays  an  important  role  in  this  period.  When 
the  battery  is  overcharged,  curves  become  much  steeper.  Heat 
generated  during  overcharging  takes  a  great  part  of  total  heat  gen¬ 
eration. 

The  curves  in  each  figure  are  quite  similar,  so  other  curves  can 
be  expressed  by  curve  of  SOC  0%  through  linear  transform: 

tsoc  =  to%  -  ^52 SOC,  SOC  =  30,  50,  70,  90,  100%, 
rate 

rate  =  1,3,  5  (1) 

where  ts oc  denotes  charging  time  of  the  battery  in  certain  SOC;  t0% 
denotes  charging  time  of  the  battery  in  SOC  0%.  Therefore,  test  result 
of  SOC  0%  is  chosen  for  curve  fitting  to  get  a  function  that  describes 


Fig.  1.  Experiment  device  system.  Each  part  of  the  device  in  which  the  tested  battery 
is  fixed  is  labeled. 
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Fig.  2.  Heat  generation  rate  of  1C  rate  charge.  Different  curves  from  the  right  to  the 
left  denote  SOC  from  0  to  100%. 


the  charge  process.  Others  can  be  obtained  by  linear  transform. 

Heat  generation  becomes  more  seriously  when  charge  rate 
increases.  One  reason  is  that  the  amount  of  heat  generated  by  side 
reaction  in  overcharging  period  remains  nearly  the  same  because 
of  the  same  starting  and  finishing  SOC  while  time  is  much  shorter. 
Another  reason  is  heat  generated  by  electrical  resistance  becomes 
considerable  when  the  current  goes  up,  which  is  known  as  P = I2R  (P 
denotes  power,  /  denotes  current  and  R  denotes  internal  resistance). 

Each  curve  is  divided  into  two  sections  to  help  us  understand 
what  happens  through  the  process: 

( 1 )  Normal  charging  section,  heat  generation  rate  is  slowly  increas¬ 
ing.  According  to  Refs.  [8-10],  the  increase  is  caused  by  reaction 
heat  and  joule  heat.  In  some  cases,  heat  generation  rate  could 
be  assumed  constant,  because  the  slope  is  small. 

( 2 )  Overcharge  section,  recombination  whereby  the  generated  oxy¬ 
gen  and  hydrogen  return  to  water  contributes  to  Q  [10]: 

2H+1o2^H20,  -285.8  kj  mol-1  (2) 

In  1C  and  3C  rate  charge,  heat  generation  rate  increase  asymp¬ 
totically  with  charging  time.  In  5C  rate  charge,  charge  time  was 
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Fig.  3.  Heat  generation  rate  of  3C  rate  charge.  Different  curves  from  the  right  to  the 
left  denote  SOC  from  0  to  100%. 
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Fig.  4.  Heat  generation  rate  of  5C  rate  charge.  Different  curves  from  the  right  to  the 
left  denote  SOC  from  0  to  100%. 


shorter  and  Q  seems  like  a  line. 
Expressions  of  Qare 


The  two  equations  are  similar  with  theoretical  equation.  Therefore, 
it  is  reasonable  to  believe  that  the  overcharge  reaction  of  1 C  and  3 C 
rate  charge  is  a  first  order  reaction  or  pseudo  first  order  reaction. 

Reaction  inside  the  battery  during  5C  rate  charge  seems  to  be 
a  zeroth  order  reaction.  It  may  be  caused  by  heat  conductance 
and  diffusion  of  reactants  inside  the  battery,  diffusion  is  the  rate¬ 
determining  step  under  high  rate.  And  the  heat  conductivity  limited 
much  heat  inside  the  battery.  Therefore,  heat  generation  rate  may 
not  help  determine  the  reaction  order  of  high  rate  charge  because 
of  the  heat  conductivity. 

Quantities  of  heat  generated  in  charge  period  are  also  measured 
by  microcalorimeter,  quantities  of  heat  of  1 C,  3  C,  and  5C  rate  charge 
are  18644.49,  24183.74,  and  27176.76 J. 

3.2.  Thermal  models 


Two-dimensional  model  is  developed  to  analyze  the  thermal 
behavior.  The  battery  is  assumed  to  be  axial  symmetric,  so  the  heat 
transport  equation  is 


V-f- 
r  dr 


(12) 


IC  :  T(r,  z,  0)  =  T0 


(13) 


For  1C  : 


f  0.74555  +  0.000158t  ,  t  <  3303 

^  -  1  11.50979  -  54303.77062  x  0.99743(  ,  3303  <  t  <  5400 

BC: 

9T 

~kr  dF 

-  h(T  r=R-T°o)> 

r=R 

z=0,H 

(3) 

=  ^TUh~T^ 

(14) 

For3C  : 


f  1 .93313  +  0.00253t  ,  t  <  1050 

\  51.09361  -  329.78716  x  0.9982f  ,  1050  <  t  <  1800  W 


For5C :  Q  = 


96826 +  0.01303t 


,  t  <  653 


-75.54081  +  0.12909t  ,  653  <  t  <  1080 


(5) 


According  to  reaction  kinetics,  reaction  rate  differential  equa¬ 
tion  of  first  order  reaction  is 

~di=kc  (6) 

c  denotes  concentration  of  reactant. 

For  the  reaction  we  presume  reaction  heat  as  AH  (kj  mol-1 ): 


d  d(AHc)  d(AHc) 
dt  dt  “  <  dt 
That  is 


and  integral  form  is 
Q  =  A  -Ae~kt 


(7) 

(8) 

(9) 


where  A  is  a  constant. 

As  for  side  reaction  mentioned  above,  expressions  of  Qfor  1C 
and  3C  rate  charge  in  overcharge  section  can  be  written  as 


For  1C  :  Q  =  11.50979  -  54303.77062  x  0.99743f 


=  11.50979  -  11.50979  x  e-0  002573308^-3287)  (10) 


For3C  :  Q  =  51.09361  -  329.78716  x  0.9982f 


(11) 


Cp  is  average  heat  capacity,  which  is  measured  by  microcalorimeter 
before  charging  and  assumed  constant;  kr  (kz)  is  thermal  conduc- 


tivity  of  the  battery  in  r-direction  (z-direction);  T0  denotes  initial 
temperature  of  the  battery  before  the  experiment;  Q  is  heat  gen¬ 
eration  rate,  which  is  expressed  by  (3)— (5);  h  denotes  heat  transfer 
coefficient.  Values  are  shown  in  Table  2. 

Too  is  ambient  temperature,  it  was  increasing  during  charging 
process,  the  value  is  estimated  by 

For  1C  :  Tx  =  24.055  +  2.71373  x  5^QQ 

(15) 

For  3C  :  Too  =  24.055  +  2.88291  x  * 

loUU 

(16) 

For  5C  :  Too  =  24.055  +  2.82556  x  ln(n 

1080 

(17) 

According  to  Newton’s  law  of  cooling: 

Q_  =  hAAT 

(18) 

A  T denotes  temperature  difference  between  battery  wall  and  ambi¬ 
ent  temperature. 

The  experiment  shows  that  AT  is  smaller  than  3°C,  so  we 
assume  it  to  be  constant  for  simplicity,  3  °C.  Heat  transfer  coefficient 
is  expressed  by 


h  = 


A  AT 


(19) 


Table  2 

Values  of  parameters  in  thermal  model 


Parameters 

Values 

Cp,  heat  capacity  (J  kg-1  °C_1 ) 

301.206 

kr,  thermal  conductivity  in  r-direction  (J  s-1  m-1  °C_1 ) 

0.74a 

kz,  thermal  conductivity  in  z-direction  (J  s-1  m- 1  °C_1 ) 

0.85a 

A,  surface  area  (m2) 

0.00769 

T0,  initial  temperature  (°C) 

24.055 

a  Ref.  [8]. 


,-0.001801621  (t-1035) 


=  51.09361  -51.09361  x  e 
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Height :  Temperature  [°C]  Max  :  28.021 


Min  :  26.766 


Fig.  5.  Temperature  distribution  by  calculation  in  tested  battery  after  1 C  rate  charge. 
Red  area  has  higher  temperature  than  blue  area.  (For  interpretation  of  the  references 
to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 

Solutions  of  three  equations  are  shown  in  Figs.  5-7,  two- 
dimensional  temperature  profile  in  the  r-direction  and  z-direction 
under  at  the  end  of  1C,  3C,  and  5C  rate  charge.  For  batteries 
charged  from  SOC  0%,  surface  temperature  of  3C  rate  charge  is 
higher  than  the  others,  a  possible  reason  is  3C  rate  charge  gen¬ 
erates  more  heat  and  there  is  enough  time  for  heat  to  get  outside 
of  the  battery.  Flowever,  for  batteries  charged  from  SOC  100%,  sur- 


Height :  Temperature  [°C]  Max  :  63.244 


Min  :  26.816 

Fig.  6.  Temperature  distribution  by  calculation  in  tested  battery  after  3C  rate  charge. 
Red  area  has  higher  temperature  than  blue  area.  (For  interpretation  of  the  references 
to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


Hight:  Temperature  [°C]  Max:  72.375 


Min:  26.741 

Fig.  7.  Temperature  distribution  by  calculation  in  tested  battery  after  5C  rate  charge. 
Red  area  has  higher  temperature  than  blue  area.  (For  interpretation  of  the  refer¬ 
ences  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the 
article.) 


face  temperature  of  1 C  rate  charge  is  higher  because  of  longer  time. 
The  battery  material  cannot  dissipate  heat  generated  immediately 
under  high  rate  charge  owing  to  low  heat  conductivity.  The  inside 
temperatures  of  3C  and  5C  rate  charge  are  much  higher  than  that 
of  1C  rate  charge,  up  to  63.24  and  72.38  °C.  The  steep  tempera¬ 
ture  gradient  may  be  caused  by  low  heat  transfer  coefficient  of 
materials  inside  the  battery,  which  is  harmful  to  the  battery  and 
would  affect  the  battery  performance.  After  1C  rate  charge  the 
inside  temperature  is  28.021  °C,  the  battery  is  under  good  working 
condition. 

Surface  temperature  rises  by  calculation  are  2.711,  2.761,  and 
2.686  °C.  Comparing  with  the  results  of  experiment,  calculation 
results  are  lower  than  actual  rises.  It  seems  that  we  have  over¬ 
estimated  the  cooling  efficiency  possibly  owing  to  estimation  of 
h.  The  temperature  rises  are  close  especially  for  1C  and  3C  rate 
charge  as  shown  in  Table  3,  so  we  believe  the  model  results  are 
reliable. 

The  temperature  distribution  is  non-uniform  because  the  poor 
conductivity  in  r-direction  limits  the  heat  transfer  during  charging 
process.  It  is  difficult  to  greatly  improve  the  heat  conductivity  of  the 
battery  because  it  is  related  to  materials  inside  the  battery  includ¬ 
ing  electrodes,  separators  and  so  on.  Therefore,  high  rate  charge 


Table  3 

Comparison  of  experiment  and  calculation 


Rate 

Temperature 

SOC 

0% 

30% 

50% 

70% 

90% 

100% 

1C 

Experiment 

Calculation 

2.714 

2.711 

2.656 

2.656 

2.436 

2.436 

2.315 

2.313 

2.289 

2.286 

2.299 

2.299 

3C 

Experiment 

Calculation 

2.883 

2.761 

2.652 

2.629 

2.460 

2.438 

2.309 

2.287 

2.093 

2.061 

1.907 

1.873 

5  C 

Experiment 

Calculation 

2.826 

2.686 

2.675 

2.642 

2.545 

2.517 

2.152 

2.126 

1.956 

1.925 

1.844 

1.794 
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should  be  avoided  in  actual  use.  It  may  cause  some  damage  to  the 
battery. 

3.3.  Conclusions 

We  check  the  surface  temperature  of  tested  battery,  calculation 
results  are  close  to  experiment  results.  The  battery  is  overcharged 
to  150%  of  its  actual  capacity,  temperature  rise  is  about  3°C.  It  is 
much  lower  than  that  under  air-cooling  condition. 

There  is  still  a  problem  that  experiment  results  and  calculation 
results  do  not  match  very  exactly  for  5 C  rate  charge.  The  tempera¬ 
ture  rise  is  getting  lower  while  heat  generation  of  5 C  rate  charge  is 
greater  than  1C  rate  charge.  A  possible  reason  is  that  in  the  model 
1C  rate  charge  experiment  takes  more  time  so  that  heat  generated 
by  tested  battery  has  enough  to  time  get  out,  and  the  temperature 
falls  down.  As  for  5C  rate  charge,  time  is  much  shorter  and  heat 
remained  inside. 
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